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Abstract 

We calculate neutrino and anti-neutrino scattering off electrons 
and nucleons in supernovae using a detailed Monte Carlo transport 
code incorporating realistic equations of state. The goal is to deter- 
mine whether particles in a neutron star core, partially spin- aligned 
by the local magnetic field, can give rise to asymmetric neutrino emis- 
sion via the standard parity breaking weak force. We conclude that 
electron scattering in a very high magnetic field does indeed give a 
net asymmetry, but that nucleon scattering, if present, removes the 
asymmetry. 
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PSR 2224+65 has a measured [0] transverse speed > 800 km/s, indi- 
cating that it has enough kinetic energy to leave the gravitational potential 
well of our galaxy, which requires a speed of 450 km/s 0. It now appears 
that many pulsars have abnormally high proper speeds Since neutrino 
emission carries off 99% of the binding energy of the magnetized neutron 
star/pulsar formed in supernova explosions, it is particularly attractive to 
consider models of asymmetric neutrino core emission as the source of the 
pulsar kick, since then only a small fraction of 0(1%) of that energy is re- 
quired to be asymmetrically radiated. In this regard, several authors have 
proposed mechanisms to accomplish this. Bisnovatyi-Kogan [|J considers 
toroidal magnetic field configurations in a differentially rotating star. In the 
absence of dissipative processes, the neutron star returns to the state of rigid 
rotation by loosening the induced higher moment fields via magnetorotational 
explosion. However, it is not very clear how this gives rise to asymmetric 
neutrino luminosities. Kusenko and Segre || , using the approximation of the 
last scattering surface method, propose that neutrino oscillations induced by 
non-zero neutrino masses alter the shape of the neutrinosphere, and derive a 
necessary tau neutrino mass of ~ 100 eV. Grasso et al || propose neutrino 
asymmetries without neutrino mass requirements, based on weak violation of 
universality. We want to study the possibility of explaining the effect within 
the Standard Model, without relying on last scattering assumptions. Instead, 
we simulate the complete process of neutrino birth and transport within the 
core. We find that the last scattering surface assumption does not reproduce 
the neutrino transport properties. While this paper was being prepared, Lai 
and Qian have considered the possibility that spin aligned nucleons could 
lead to asymmetric neutrino emission. We, in fact, find just the opposite: 
nucleons are parasitic and it is the electrons which have the potential to lead 
to asymmetric emission. 

The spin-aligned fraction, r] = [N(+) - N(-)]/[N(+) + N(-)\, for non- 
relativistic non-degenerate particles (i. e. nucleons in proto - neutron star 
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cores) with magnetic moment /im in a magnetic field B at temperature T is 

r/ = tanh^^, (1) 

while for relativistic ideal degenerate particles (i. e. electrons in neutron star 
cores) (with T << Ep, Ep = Fermi energy), 

_ {Ep + (imB) 3 — (E F — fiAjB) 3 
V ~ (Ep + m Bf + (E F -^ M Bf U 

A cooling proto-neutron star core at T = 5 MeV, B = 10 14 Gauss, with baryon 
density of 5 x 10 38 baryons/cm 3 , 3% of it as protons, has ?7neutron = 1-2 x 10~ 4 
and r/eiectron = .011. Thus nucleons are hardly aligned at all, and electrons, 
due to the relativistic degeneracy, have also a rather small spin alignment. 
Only for gigantic fields B ~ 2000 GT will the ideal gas electrons be highly 
spin aligned in neutron stars. It is possible, however, that electrons do not 
form an ideal relativistic degenerate gas within neutron stars, but instead be- 
come nearly a perfect ferromagnetic spin aligned fluid, with a nearest neigh- 
bor Hamiltonian ~ YUj &i ' Sjhij w hh h it j < 0. This may happen since 
parallel spin alignments will tend to reduce the Coulomb repulsion between 
electrons, a major energy consideration at these high number densities. If 
such an electron ferromagnetic fluid were present, then electron spin align- 
ment would be nearly 100% even for 10 14 Gauss. In our simulation, we just 
use an ideal electron gas, disregarding possible ferromagnetic effects. 

There has been some recent discussion in the literature indicating that no 
asymmetry can arise from neutrino multiple scattering if the neutrinos are in 
thermal equilibrium ||. We fully agree with this statement, although it is not 
the case we are considering here. The neutrinos that produce the asymmetry 
do not reach thermal equilibrium, due to a non-negligible neutrino mean free 
path comparable to the temperature gradient scale. In fact, too little or 
too many collisions, leads to zero neutrino asymmetry, but we find that a 
moderate number of collisions per neutrino of ~ 4000 leads to the effects 
reported here. 
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n-v 


p-v 




e-v e 




c v 


-1/2 


(l-4sin 2 6 w )/2 


L e 


2 sin 2 9 W + 1 


2 sin^ 9 W - 1 


C a 


-9A/2 


9a/2 


Re 


2 sin 2 9w 


2 sin 2 6^ 



Table 1: The couplings c v and c a for nucleon scattering (Eq. 7), and L e and i? e 
for electron scattering (Eqs. 4-5). For nucleon scattering with antineutrinos, 
c a changes sign; for electron scattering with neutrinos or antineutrinos the 
couplings are the same. We use sin 2 9w = 0.23 and gA = 1.26. 

The probability P(9 4 , (f) 4 )d9 4 d(f) 4 that a neutrino or anti- neutrino will scat- 
ter into the solid angle dQ 4 = sin 9 4 d9 4 d(p 4 is 

\M(pi,p 2 ,p 3 ,p 4 )\ 2 sin9 4 



P(0 4 , <Pi)d9 A 



-d9 A 



(3) 



/ \M(p 1 ,p 2 ,p 3 ,p 4 )\ 2 dtt 4 
where M(pi,p 2 ,p 3 ,p 4 ) is the scattering amplitude for a (anti) neutrino with 
incoming momentum p 2 and outgoing momentum p 4 , off a target with initial 
and final momenta p\ and p 3 respectively. For scattering of neutrinos off 
electrons, the spin-dependent matrix element squared is 

\M[e x v 2 -> e 3 z/ 4 )| 2 = 16G 2 F {Ll[(p 4 • p 3 ){p 2 • Pi) - m e (p 4 - p 3 )(p 2 ■ s x )] 

+i? 2 [(p 4 -Pi)(P3 -ftj) +m e (p 3 -p 2 )(p 4 ■ s^) (4) 
-L e R e [m 2 e (p 4 ■ p 2 ) + m e (p 2 ■ p 1 )(p 4 ■ s x ) - m e (p 4 -p 1 )(p 2 ■ s x )\} , 

where the weak couplings L e and _R e are shown in Table 1, and si is the 
polarization four-vector of the initial electron, which satisfies p\ ■ s± — and 
si = —1. For antineutrino scattering, the respective expression is the same 
as the one above after exchanging p 2 <->• p 4 . 

For nucleon scattering, Ref. [0] has given the scattering probability 



74, <p 4 



4 Uif> 4 



/do 



d0 



4 U^4, 



(5) 



where the cross section for neutrino-nucleon scattering is: 

r 2 P 2 

' ^-{c 2 +3c 2 + (c 2 -c 2 ) cos6»+2^c a [(c,;-c a ) cos^ out +(c„+c a ) cos# iri ]} 



da GlEl 



dVt 4 



4tt 2 



(6) 

Here r\ is the polarization factor of Eq. (1), c v and c a are given in Table 1, 9 
is the scattering angle, and 9i n and 9 out are the angles between the nucleon 
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Star Parameters 


B 


0.1 -> 5000 [GTesla] 


T 


1 -> 100 [MeV] 


n B 


5 x 10 37 — > 1.16 x 10 39 [baryon/cm 3 ] 


fe 


0.001^0.5 



Table 2: The physical parameters involved in the asymmetry calculation. 

spin and the incoming or outgoing neutrino momentum, respectively. In 
Table 2 we give the star parameter ranges used in our simulations, including 
the assumed large magnetic fields associated with magnestars. 

We now add the astrophysics. The density of matter in the interior of 
neutron stars increases from 10 25 baryon/cm 3 in the surface layer to ~ 10 39 
baryon/cm 3 near the center. However, more than 95% of the mass of a typ- 
ical neutron star is compressed to baryon densities exceeding 10 38 cm -3 [10 
and, to an excellent approximation, the supernova core can be described as 
hadronic matter at constant density. The different equations of state (EOS) 
give rise to different fractions of electrons present, f e = n e /riB, where n e and 
tib are the electron and baryon densities respectively. Thus we describe the 
EOS landscape by two parameters: n B and f e . We set the mass of the core to 
1.8 Mq. The thermodynamic parameters are the temperature T and the inte- 
rior magnetic field. We use constant interior magnetic fields. With respect to 
neutrino production, we only consider neutrino-anti-neutrino thermal pairs, 
as 90% of the neutrino emission is of thermal origin. 

A Monte Carlo neutrino transport code for neutrino scattering off po- 
larized targets (electrons only or both electrons and nucleons) was written 
to determine if the parity violation in the weak force can lead to asymmet- 
ric neutrino emission from the surface, and thus possibly explain the large 
proper motion of pulsars. Because of Fermi degeneracy, essentially only those 
electrons near the Fermi level participate in scattering. Since the electrons 
do not interact via the hadronic force (and the star is electrically neutral) 
they might be described as an ideal relativistic degenerate Fermi gas with 
Fermi momentum pp = 6.10464 x 10 _11 (/ e n^) 1 / 3 MeV/c, with rig in 1/cm 3 
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(but see above qualifications concerning possible electron ferromagnetism). 
An experimental constraint on the EOS is that the radius of the neutron star 



must be smaller than 14 - 16 km 11 



The simulation goes as follows: a random flavor v-v pair is born at a 
random location within the core, with a spherically random direction and 
each particle with an energy equal to 3T, where T is the temperature at 
that position in the star. A discrete temperature gradient was imposed. 
For this purpose, four concentric shells of inner and outer radiae \R a tw and 
^^--Rstar respectively (n = 0,1,2,3) were considered. The temperature in 
each shell was set to ^pT.,, where T s is a temperature scale parameter. Thus 
for birth location in the nth shell the neutrino energy was E U)P = (6 — n)T s . 
For a uniform temperature in a constant density core, the v-v pair emission 
probability per unit volume is also uniform; the production point thus has 
a radius r ~ R sta r X £^ 3 , where £ G [0,1] is a uniformly distributed random 
variable. Instead, with negative temperature gradients (hotter center region), 
the emission is also higher at the center, skewing the distribution towards 
r ~ Rstar x £• We used here r = 0.4i? star x £, with the phenomenological 
factor 0.4 to take into account v-v production only within the interior core. 
After birth, the location of the next scattering event is statistically generated 
from the distribution given by the mean free path I = l/(Y!,i n i a i) i n i is the 
i-type target density -electrons or nucleons- and cr, the cross section) |TJ . 
If this is outside the radius, the neutrino energy and momentum components 
are accrued as ejected flux. If the scattering position is still within the core, 
a randomly directed electron (in the case of electron collisions) from the 
Fermi surface is picked, with its spin direction parallel or anti-parallel to the 
magnetic field, determined from the spin probability distribution function. 
The probability of parallel/antiparallel spin is defined in terms of r\ from 
Eq. (2) as P± — (1 ± i])/2. The direction of the outgoing neutrino is then 
generated according to the distribution given by the differential cross section 
(Eq. 3) and a Pauli check is made. Then the outgoing neutrino seeks the 
next scattering event. This is continued until the neutrino leaves the star. 
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Finally, the history of the neutrino's birth partner is followed. This is done 
for N neutrinos. 



The asymmetry is defined as the accrued ejected momentum over the 
accrued ejected energy: 



Asymmetry = =^ , (7) 



where z is the direction of the magnetic field ||13|| , We verified that the other 
momentum components, J2p x and J2p y , fall within the statistical fluctua- 
tions, while J2p z stands out as the signal. The variance of the statistical 
fluctuations increases as ~ v~N, so the asymmetrical statistical 'background' 
decreases ~ 1/yN. To obtain a 0(1%) positive signal with a ~ .1% error 
thus requires N ~ 10 6 . In the actual runs, one million pairs were followed, 
giving iV = 2 x 10 6 . The 'background' statistical asymmetry was calculated 
by performing a run with random spin orientation at each scattering location, 
which, for the samples of a million pairs gave a typical magnitude of 0.05%. 
For runs including both electron and nucleon scattering, the probability of 
electron scattering is -P e iectron = S e i ec t r on/ (^electron + ^proton + S neu tron)) where 
Sj = rii (particle density) xo~j (scattering cross section). Similar probabilities 
exist for scattering off neutrons and protons. 

The results we present in our tables and figures correspond to a neutron 
star with 1.8 solar masses, radius 8.3 km and electron/nucleon ratio f e = 
0.12. This size and structure is representative of most realistic equations of 
state which incorporate nuclear repulsion at short distances. Typically, the 
heaviest neutron star for such equations of state is over two solar masses (see 
e.g. R. Bowers et al. in Ref [10]). 

Several critical details arose. Firstly, it is a highly non-trivial problem in 
computational physics to draw scattering samples from Eqs. (3-4). We use 
the 'rejection method' fl~4" |, generalized to two random variables. Secondly, 



due to the billions of random numbers used, we discovered that a popular 
random number generator, Knuth's algorithm flT5|] , performed unsatisfacto- 
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T 5 (MeV) B (GT) target particles 


Asymmetry 


10 200 electrons 
10 4000 electrons 
10 4000 nucleons & electrons 


(0.18 ±0.05)% 
(0.47 ±0.05)% 
(0.04 ±0.05)% 



Table 3: Asymmetry in neutrino flux for different magnetic field and scatter- 
ing targets, for a neutron star of 1.8 M , radius 8.3 km, and electron/nucleon 
fraction 0.12 

rily. We finally used the 'gold-plated' random number generator ran2, of 
Ref. fl6f , which performed without problems. 



In Fig. 1 we display the results of a positive signal, corresponding to the 
case of extremely high magnetic field (~ 10 3 GT) and scattering with the 
electron gas only. Subtle physics effects are at play here. First, looking at 
Eqs. (3-4), we see that the spin terms are multiplied by m e , the electron 
mass; thus, it would seem (wrongly) that relativistic electrons have little 
parity violation. Secondly, from Fig. 1, the asymmetry increases with tem- 
perature, reaching 0(1%) at ~ 17 MeV. These two observations are related: 
the magnitude of the polarization four-vector for electrons traveling paral- 
lel or anti-parallel to the magnetic field is proportional to the relativistic 7 
factor and so for such electrons m e s ~ E; as the temperature increases, the 
Pauli window increases allowing a greater energy range of allowed transitions. 
Thus the asymmetry is expected to be an increasing function of temperature. 

Now, if we include neutrino scattering with nucleons as well as with elec- 
trons, the results change dramatically. Elastic scattering of neutrinos off 
hadrons occurs via the neutral current. The small magnetic moment of the 
nucleons gives very little polarization, so we would expect that allowance 
of nucleon scattering would tend to remove the neutrino emission asymme- 
try. Moreover, the nucleon density is much larger than the electron density, 
so that nucleon scattering and its effects should dominate. In Table 3, we 
present the results of the nucleon scattering inclusion, showing that this is 
indeed the case: even at extremely high magnetic fields, the inclusion of 
nucleon scattering washes out the asymmetry, as shown in the last row of 
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Table 3. Here we should recall that the statistical fluctuation of the simu- 
lations are about 0.05% in the asymmetry, which is already larger than the 
signal obtained. 

Our numerical results support the hypothesis that neutrino scattering 
off spin-polarized electrons is a viable mechanism for achieving high pulsar 
proper velocities in the case where extremely high magnetic fields exist in 
the core, and whenever scattering with nucleons can be neglected, without 
resorting to any non-standard physics or neutrino oscillation mechanism, but 
if nucleon scattering is present, the asymmetry disappears. Now, since the 
nucleon density is necessarily larger than the electron density, to neglect 
nucleon scattering in favor of electron scattering does not seem plausible, 
unless there is an additional mechanism, like superfluid phases in the core, 
that could block nucleon scattering. 
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Figure 1: Asymmetry vs. Temperature, for a star with baryon density % = 
8 x 10 38 1/cm 3 , radius -R sta r = 8.3 km, electron fraction f e = 0.12, and 
magnetic field B = 4000 GTesla. 
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